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Abstract 

We compare the classical and quantum mechanical position-space probability 
densities for a particle in an asymmetric infinite well. In an idealized system with 
a discontinuous step in the middle of the well, the classical and quantum proba- 
bility distributions agree fairly well, even for relatively small quantum numbers, 
except for anomalous cases which are due to the unphysical nature of the poten- 
tial. We are able to derive upper and lower bounds on the differences between 
the quantum and classical results. We also qualitatively discuss the momentum- 
space probability densities for this system using intuitive ideas about how much 
time a classical particle spends in various parts of the well. This system provides 
an excellent example of a non-trivial, but tractable, quantum mechanical bound 
state problem where the correlations between the amplitude and curvature of 
quantum mechanical wavefunctions can be easily compared to classical intuition 
about particle motion, with quantitative success, but also warning of possible 
surprises in non-physical limiting cases. 
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1. Introduction 

The calculation of position-space energy eigenstates in quantum mechanical bound 
state systems is only the first step in student understanding of the underlying physics. 
Very often, a comparison of the resulting quantum probability densities, versus 
X, with classical probability distributions provides much needed insight into the deep 
and pedagogically useful connections between the two approaches to mechanics. The 
most familiar problem of one- dimensional quantum theory, namely the infinite well, is, 
for example, occasionally discussed in this manner [1, 2]. The local averaging required 
to show that the oscillatory quantum probability distribution approaches the well- 
known flat classical value is easily implemented in this system, either visually or more 
analytically, to illustrate in detail the correspondence principle limit. 

A far more familiar example and image in many modern physics and quantum me- 
chanics texts is the comparison of the classical and quantum results for the harmonic 
oscillator eigenstates for increasingly large quantum number, appearing in both older 
classic texts [3, 4] and as well as in very recent ones [1, 2], [5] - [10] at a wide vari- 
ety of levels. In this case, both the amplitude and the curvature ('wiggliness') of the 
wavefunction vary in a non-trivial and highly correlated way which can be understood 
using classical connections. A number of texts [1, 2], [6], [11] - [13] now emphasize the 
intuitive ideas behind the form of wavefunctions as dictated by the shape of the poten- 
tial energy function and the value of the quantized energy eigenvalue, with examples 
including less familiar systems such as linear potentials [2, 8, 13, 14] and asymmetric 
infinite wells [2, 12, 13]. Similar comparisons are also possible for two-dimensional sys- 
tems (such as the circular infinite well [15]) and an excellent discussion of the classical 
limit of the quantum solutions for the hydrogenic radial probability distributions has 
appeared in the pages of this journal [16]. 
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Most such presentations can easily give the impression that the approach from the 
case of low- lying quantum states states (small n) to the classical limit (n >> 1) is 
smooth and uninteresting and in this note we wish to examine a simple system which 
exhibits some unexpected properties. The potential we consider is an asymmetric 
infinite well defined via 



V(x) = < 



' oo for x < —a 
for —a < X < 



Vo for < X < +b ^"^^ 



^ oo for +6 < a; 

so that it is an infinite well of width {a+b) , but with the right side at a higher (constant) 
value of potential. (In most of our numerical calculations, we will actually use values 
of a = 6, but we will analyze the problem in some generality, at least initially.) Such 
a potential is shown in Fig. 1 (with specific values of a, 6, and Vq for later use.) This 
potential is useful for the discussion of semiclassical limits for several reasons: 

(i) It is a simple example where there is a non-trivial but easy-to-visualize variation 
of potential energy, and hence speed, between the classical turning points, 

(ii) the classical concepts of how much time a particle spends in each side of the 
well are intuitively obvious so that the classical probability distributions are 
straightforwardly obtained, 

(iii) the quantum solutions can be obtained in closed form, and the resulting energy 
eigenvalue conditions implemented numerically rather simply, and 

(iv) the Fourier transform of the position space solutions can be readily obtained to 
discuss the same qualitative and intuitive ideas in momentum space. 

It is not surprising, perhaps, that several textbooks [7], [12], [13] use this system as a 
qualitative example of intuitive wavefunction analysis. 



2. Quantum and Classical Solutions 



To obtain the quantum solutions for this problem, we assume, at least initially, 
that E > Vq (such as for those states labeled 5 and above in Fig. 1), and we solve the 
time-independent Schrodinger equation in each side of the well, obtaining solutions of 
the form 

. , \ _ j ^ sin[A;(a: + a)] for —a < a; < . ^ 

~ \ B sm[q{x - b)] ior < x < +b ^ ' 

where 



and (3, 

These solutions already satisfy the boundary conditions at the two infinite walls. In- 
sisting on the continuity of ijj and ijj' at the origin, the energy eigenvalue condition is 
given by 

k cos{ka) sm{qb) + q cos{qb) sm{ka) — (4) 

which is easily solved graphically and/or numerically. In the case where < Vq, we 
can let 

/2m(E - Vo) . /2m(Vo - E) ._ , . 

^ = V — ^^V — — 

and use the same solutions and eigenvalue condition, but with the substitutions 

sin(g6) —>■ sm{iqb) = i smh.{qb) 

cos{qb) — > cos{iqb) — cosh(g6) (6) 

and many multi-purpose mathematical packages such as Mathematical can easily ac- 
commodate this change automatically. The energy eigenvalues can be generated for 
any given set of numerical parameters and the corresponding wavef unctions in (2) 
obtained and normalized. 
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For purposes of comparison with purely classical results, we note that the classical 
probability distribution for the particle when E < Vq (where it would be restricted to 
bounce back and forth between the walls at —a and 0) is 

_ , . f 1/a for —a < x < Q 
Pcl{x) = I q for < o; < +6 ■ 

This imphes that the probability of finding the particle in the left side of the well in 
this case is 

Prob[-a <x<0] = Pf^'^ = 1 for £; < Vq. (8) 

For the case when £' > Vq, the situation is more interesting and we can approach it 
by first calculating the time spent (classically) by the particle in the left (L) and right 
(i?) sides of the well, namely 



Tj, = 2a/vL = 2a 



m 
2E 



These combine to give the (classical) probabihty of finding the particle in the left side 
as 

(CL) _ _ q/VE 



Prob[-a < a; < 0] = P^^^' = 



Tl + Tr aj^fE + hj^E - T/o 
a^JE - Vq 

~ ay/E - Vq + h^/E 

a 



(10) 



a + h^E/{E - Vq) 

ior E > Vq. (The final form is motivated by a discussion below.) In a similar way we 
find that 

Prob[0 < X < +6] ^ Pf^) = fw^^i^fw - (11) 

a^/E — Vo + by E 
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In the limit where E » Vq (and the higher potential on one side has little effect) we 
have the purely geometric results 

piCL)^a Pf^)^-^. (12) 

^ a+b ^ a+b ^ ' 

Combining these results, we can evaluate the classical probability distribution for 
the general E > Vq case, namely 

I ^E- Vo /{aVE -Vo + b^/E) for -a < x < , . 

^CL[x) I ^^^^^^^—^^ ^ for < X < +6 ^ ^ 

which is properly normalized since 

PcLix)dx = l. (14) 



In order to compare these classical quantities to the quantum results, we will eventually 
require the normalized, position-space energy eigenstates to evaluate 

P^Q^^ ^ r dx and P^^""^ = \iP{x)\^ dx. (15) 

J-a Jo 

As an example of the possible solutions in such a system, we first choose a standard 
set of parameters, namely 

n^2m^l and a = 6 = 3 , Vo ^ 20 (16) 

and examine the results in detail, indicating below how general they are: specifically, 
we choose equal values of a and b to facilitate comparison between the quantum and 
classical probabilities of finding the particle in either side. With these parameters, the 
lowest 9 energy eigenvalues are given by 

K = 0.95 , 3.78 , 8.44 , 14.78 , 20.84 , 22.34 , 24.94 , 29.24 , 33.30 (17) 

for n — 1, ...,9 and these are the values shown in Fig. 1 (as horizontal dashed lines) 
for the standard parameter set. (We note that the first four values, corresponding to 



states mostly localized to the left side, have values which are appropriately smaller 
than the first four states for an infinite well of width a, namely = ^^7r^n^/2ma^; 
with this parameter set, for example, we find that E^/ Ri 0.85 < 1 since the 
wavefunctions can penetrate into the classically disallowed region, thereby reducing 
their 'wiggliness' and kinetic energy.) We next plot, in Fig. 2, the normalized position- 
space probability densities, versus x, for these states, using the appropriately 
normalized eigenf unctions. The vertical dashed line indicates the center of the well, 
while the sohd horizontal hues (one only for the E < Vq cases for 1 — 4, and two for 
the E > Vq cases for 5 — 9) are the classical probabihty distributions in (7) and (13) 
respectively. 

We first note the familiar general pattern of the addition of one 'effective half 
wavelength' or additional node as the quantum number increases, typical of all such 
bound state systems. For the 1 — 4 states, which are all below the Vq threshold, the 
similarity to the corresponding infinite well results is clear, with the additional feature 
of the increased tunneling into the classically disallowed region. The first state above 
threshold, n = 5, is consistent with the usual intuitive rules about the behavior of 
wavefunctions in different potential regions. The wavefunction is wigglier (less wiggly) 
but with lower (higher) amplitude in the left (right) side where the classical kinetic 
energy and speed is larger (smaller) and the classical particle would expected to spend 
less (more) time. The n = 5 state is even roughly consistent with the classical result 
(13), as far as such a low- lying state can be. Similar expected patterns (and qualitative 
agreement) are evident in the n — 7 and n = 8 states. 

However, we note that the n — 6 wavefunction has an unexpected form, at least in 
terms of the relative values of the amphtudes on the left and right sides. In this case (as 
well as that for n = 9 where a similar effect is seen to occur) the wavefunctions being 
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matched at the (discontinuous) x = boundary are connected (smoothly, of course) at 
what happens to be very close to an antinode. Since the amplitude in each side of the 
well is constant (since the potential is piecewise constant) , if the wavefunctions match 
at an anti-node, the amplitudes must match {A = B in (2)) everywhere in the well. 
Thus, while the 'wiggliness' of the wavefunction agrees with naive expectations based 
on considerations of classical speed and/or kinetic energy, the amplitudes and resulting 
probability densities are somewhat anomalous in that the particle will be found much 
more frequently (about half of the time in this case) in the left side of the well than 
expected from purely classical arguments. (Presumably, a student asked to sketch a 
possible solution in this potential well, as in Refs. [7, 12], and who provided something 
hke this n — 6 case would not have received full credit!) 

To discuss these effects more quantitatively, we calculate the quantum probability 
of finding the particle in the left-half of the well (15) for a number of states in this 
model system for comparison with the classical results given by (8) and (10). We plot, 
in Fig. 3, the quantum values (diamonds), pj^^^\ for all states with E < 100 (up to 
n = 18 in this case) for comparison to the solid curves representing the classical result. 
The horizontal dashed line indicates the high energy limit of Pl = 0.5 (when a = b) 
where the presence of a small potential 'bump' at the bottom of a very deep well (when 
E » Vq) will have little effect. We first note that for the states below threshold, 
the quantum probabilities, p^^\ decrease from near the purely classical value of 
Pl — 1 due to the increasing amount of quantum tunneling evident in Fig. 1. Above 
threshold, the quantum probabilities seem to track the classical prediction (10), except 
most dramatically for the anomalous cases like n = 6, 9. In these instances, pj^^^ is 
much larger than expected classically, but never exceeds Pl — 0.5 which would truly 
be unphysical as it would imply (classically) that the particle is moving slower in the 
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left side of the well and therefore spending more time there. A similar, but somewhat 
less dramatic, effect is also evident for those states with quantum probabilities which 
are significantly less than the classical predictions, such as the n — S state in Fig. 3. 
In all cases we have studied where Pl is most obviously much smaller than classical 
expectations, we have found that the wavefunction matching occurs very close to a 
node. The two extreme cases of anomalously large (small) values of pj^'^^ have thus 
been found to be connected to wavefunction matching near antinodes (nodes) at the 
discontinuous boundary. 

Using this observation, we have been able to derive upper and lower bounds for 
pj^^^ which we have found to be satisfied in all of our numerical studies. For example, 
in the case of matching at an antinode, the boundary condition on the left/right am- 
plitudes (from the continuity of ip) is A = ±B and the resulting quantum probabilities 
can be evaluated to find that 

- ^ ^^^^ 

which is also the high energy {E » Vq) limit. For the case of wavefunction matching at 
nodes, the appropriate condition comes from the continuity of namely kA — ±qB, 
which then gives 

p{QM,min) _ <^ /,qv 

^ a + b[E/{E-Vo)y ^ ^ 

(These results are considerably more simple in form than the general expression for 
piQM) gjj^^^g integrating over integral numbers of half- and quarter-wavelengths 

in each side of the well in these special cases.) These two forms can be compared to 
the purely classical result (10) and we note that 

(node case) ^--H < pf ^) < p^max) (antinode case) (20) 
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since ^jEjiE — Vq) > 1. To illustrate these bounds, we have plotted p^*?^'"*"^) (hori- 
zontal dashed line) and pj^^'"^^"-^ (dotted curve) in Fig. 3 where they clearly bracket 
the classical result. We find that the results presented there are very typical of all 
the cases we have studied, namely that the quantum results for p^^'^ rigorously fall 
between the upper and lower bounds, often coming very close to saturating them. 

3. Quantum Results in a Smoothed Asymmetric Infinite Well 

An obvious question about these effects (especially the anomalously large values 
of pj^^^^ is whether they are artifacts of the specific set of model parameters (16) 
used here or a more general phenomenon related to the structure of the potential well 
itself. We have repeated our analysis for a large number of different values of a, 6, and 
Vq (including many cases where a ^ b) and have almost always found that in the first 
10-15 eigenstates above the Vq threshold that there are 1-2 solutions which exhibit 
the wavefunction matching near antinodes which lead to anomalously large values of 
piQM)^ compared to the classical result (10), almost saturating the upper bound (18), 
as well as 1-2 states which come close to saturating the lower bound (19). We think 
that it is easy to argue that these effects are a result of the discontinuous (and hence 
unphysical) nature of the potential step ai x — 0. In order to test this, we have solved 
the Schrodinger equation for a 'smoothed' version of this asymmetric well given by 



oo for a; < —a 

Vo/{l + e-^/-^) for -a<x<+b . (21) 
oo for +b < X 



This version still has impenetrable walls at x — —a, +b, but gives a smoother transition 
between V{—a) ~ and V{+b) ~ Vq depending on the value of S: the discontinuous 
step potential is recovered in the limit that S ^ 0. This smoothed version is shown 
in Fig. 4 (middle) for 6 — 0.2 (as the dashed curve) along with the discontinuous 



{5 = 0.0) step potential. We then solve the Schrodinger equation (numerically) to 
find the allowed energy eigenvalues and the resulting normalized eigenstates for the 
smoother case. The results for the n — 6 (top) and n — 7 (bottom) states are also 
shown in Fig. 4 where the wavefunctions for the 6 = 0.2 cases arc shown (dashed 
curves) for comparison to the discontinuous case (solid curves.) For the anomalous 
n — 6 case, even the introduction of a small bit of 'smoothness' into the potential 
allows the wavefunction to accommodate to the required matching and the resulting 
quantum probability density is much closer to the classical expectation. On the other 
hand, for the n — 7 case, which was already in fairly good agreement with the classical 
result, the changes are much less dramatic. Similar smoothing functions yield similar 
results as we can see by considering a linear smoothing potential given by 

{0 for —a < X < —e 

Vo{l + x/e)/2 for-e<x<+e . (22) 
Vq for +e <x < +b 

If we expand both potentials (21) and (22) near the origin {x a; 0) we find that they 
have a similar form, namely 

For numerical purposes, therefore, we will use e — 2S for comparisons between the two 
smoothings. For example, in Fig. 4, we include the linear extrapolation (center) and 
the results for the n — 6 (top) and n — 7 (bottom) states as dotted curves (for e = 0.4). 
Not surprisingly, perhaps, the results are similar to the exponential smoothing. 

In order to quantify the improvement in agreement with the classical results, we 
show in Fig. 5 the quantum probabilities, p^^\ for both the discontinuous {5 — e — 
0.0, diamond data as before) case and a smooth exponential {5 — 0.2, starred data) 
case, compared again to the classical result (10). The agreement, in general, is much 
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better for the more physical potential gradient, except very near threshold (the n = 5 
case) where small changes in the shape of the potential might be expected to have 
large results. We note that in each case (except for n — 5) that there is improvement 
(often substantial) in the agreement with the classical result. 

At the same time, we show in Fig. 5 (top) the fractional change in the energy 
eigenvalue for the new smoothed potential compared to the discontinuous case, namely 
AE/E, and note that while the energies are different by no more than 3% in the 
smoothed potential, the changes in the corresponding values of p^^^ are much more 
substantial, up to 100%. This imphes that the introduction of the smooth (and hence 
more physical) potential step does not change the intrinsic properties of the system 
(energies, etc.) very much, but does allow for a more realistic description of the 
quantum probability densities and their approach to the classical limit. We find that 
in almost all cases where this type of anomalous behavior is encountered that the 
introduction of smoothing with a length scale (5 ~ e ~ A/4 (where A is the wavelength 
in the left side of the well) is enough to substantially improve agreement with the 
classical probability predictions while making little change in the energy eigenvalues. 

4. Momentum Space Results 

We can also develop our semi-classical intuition, as well as observing some of the 
anomalous quantum behavior, by examining the momentum-space probability densities 
for this system. The momentum-space wavef unctions, 0(p), can be obtained from the 
■0(x) in (2) using the Fourier transform via 



For the 5 — 0.0 case, where we can use the analytic results (2), the resulting normahzed 
momentum-space probability distributions are easy to generate and we show them for 




(24) 
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the same first 9 eigenvalues in Fig. 6. The vertical dashed lines correspond to the values 
of p = ±hk (for all values of E) while the dotted lines indicate values of p = ±hq (for 
states with E > Vq) where we would expect to see features corresponding to classical, 
back-and-forth motions in the left side {p = ±hk) and right sides (p = ±hq) of the well. 
For the states below threshold, the results are similar to familiar ones [2, 14] for the 
standard infinite well. For the n — 5 case, the first above threshold, we see a central 
feature, much like that for n 1, as well as two small features at larger values of |p|. 
The similarity between the Ei = 0.95, k = 0.974 central feature for n = 1 and the 
E5 — V0 — 0.84, q — 0.914 central peak for n = 5 is clear. The smaller features at larger 
values of p = ±hk are qualitatively consistent with the small amount of time spent 
in the left-side of the well for this case where q « k. The anomalously large value 
of p1^^^ ^ 0.5 for the n = 6 case, where the particle spends roughly equal amounts 
of time in each side of the well, is evident in the approximately equal magnitudes of 
the q, k 'bumps' in this case. We note that in general we expect to see slightly larger 
(smaller) features at the smaller (larger) values of p = ±?iq (p = ±hk) due to our 
intuitive probability arguments (Pr > Pl) and this pattern is also reasonably evident 
for the n = 7, 8 cases. This pattern becomes increasingly difficult to identify, however, 
as n increases because the k, q features tend to merge due to the fact that 



which goes to zero as E increases. The difficulty with unambiguously identifying them 

is already evident in the n = 9 case where the classical k, q features are obvious, but 
the interference between the k, q terms gives rise to a feature between them which is 
even larger. 




(25) 



2y/E 
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It is tempting to imagine trying to correlate, in a quantitative manner, the 'amount 
of probability' in the ±q and ±k 'peaks' with the position-space probabilities of mea- 
suring the particles in the right and left sides of the well respectively. We find, however, 
that given the large amount of interference between the k, q pieces of the Fourier trans- 
form necessary to obtain (p{p) that no such identification is possible. 

5. Conclusions 

In conclusion, we have identified a simple model quantum mechanical bound state 
system, the asymmetric infinite well, which can be easily analyzed and for which the 
comparison of the quantum and classical probability distributions finds many points of 
similarity, even exhibiting quantitative agreement for probabilities for many states. We 
have even been able to derive seemingly rigorous upper and lower bounds on how much 
the quantum results deviate from the classical probability expectations. The system 
can also be used to qualitatively examine classical intuition about momentum-space 
probabilities as well, given the relative simplicity of the position-space wavef unctions. 
However, it seems to yield, in almost every case studied, some solutions which have 
anomalous probability distributions, compared to classical expectations, due to the un- 
physical (discontinuous) nature of the potential. It also serves, therefore, as something 
of a cautionary tale about the use of idealized models. 
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Figure Captions 



Fig. 1. Example of the asymmetric infinite well potential (1) for the special case of 
a = b = 3 and Vq = 20. The nine lowest energy levels, obtained from the energy 
eigenvalue equation (4) with the parameters in (16), are shown as horizontal 
dashed hues. 

Fig. 2. Normahzed position-space probability densities, versus x, for the first 

nine energy levels in the asymmetric well shown in Fig. 1 with the parameters 
in (16). The sohd horizontal hues indicate the classical probability distributions 

(7) (for E <Vo and states 1 - 4) and (13) (for E > Vq and states 5-9). The 
vertical dashed lines indicate the center of the well at x = 0. 

Fig. 3. Classical probability of finding the particle in the left side of the well (solid line 
and curve), given by (8) (for E < Vq and states 1—4) and (10) (for E > Vo 
and states 5 — 9), as a function of E for the states in Figs. 1 and 2 with the 
parameters in (16). The data points (diamonds) are the corresponding quantum 
probabilities, pj^^^ (15)^ for the lowest 18 energy eigenvalues (all those with 
E < 100) with the standard parameter set. The horizontal dashed line at Pl — 
0.5 corresponds to equal amounts of time spent in the left and right halves of 
the well which is also the upper bound (18), Pf'^'™""^ for this case. The dotted 
curve is the corresponding lower bound (19), pjQ^'^^"\ 

Fig. 4. The middle figure shows two 'smoothed' versions of the discontinuous step (solid 
curve) at the center of the well given by (21) (S — 0.2, dashed curve) and (22) 
(e = 0.4, dotted lines). The value of e = 25 is chosen so that the two smoothings 
as in (23). The solutions corresponding to n = 6 (top) and 
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n = 7 (bottom) for the discontinuous well {5 = e = 0.0, solid curves) and the 
smooth cases with S = 0.2 (dashed curves) in (21) and e = 0.4 (dotted curves) in 
(22) are shown. The 'anomalous' n — 6 case is affected dramatically and becomes 
much more consistent with the classical predictions for the smoothed cases, while 
the more standard n — 7 case is not affected very much. 

Fig. 5. Same as Fig. 3, but for £^ > Vq only. The solid curve is the classical prediction 
for Pj^^\ and the results for the quantum probability from (15) for the discon- 
tinuous case {S = 0.0, diamonds again) and the smoothed case (21) with 6 = 0.2 
(stars) are shown. The upper (dashed) and lower (dotted) bounds on pf^^^ from 
(18) and (19) are indicated. Even a small amount of smoothing of the discon- 
tinuity improves the agreement with the expected classical result dramatically. 
The fractional change in energy t^EjE induced for each state by using 5 = 0.2 
compared to 5 = 0.0 is shown at the top. In no case is the change larger than 
3% and it decreases quickly as the energy is increased. 

Fig. 6. Normalized momentum-space probability distributions, |0(p)P versus p, for the 
lowest nine states in the asymmetric well with the standard parameter set (16). 
The vertical dashed lines show values of p = ±^A; (for all E) and the dot- 
ted lines correspond to p = -^Uq (for E > Vq) from (3) for comparison. The 
momentum-space wavefunctions are calculated by taking the Fourier transform 
of the position-space solutions in (2). 
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